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Abstract 

Since the discovery of the first exoplanet, those most adequate for life to begin 
and evolve have been sought. Due to observational bias, however, most of the discov- 
ered planets so far are gas giants, precluding their habitability. If these hot Jupiters 
are located in the habitable zone of their host stars, and if rocky moons orbit them, 
these moons may be habitable. In this work we present a model for planetary transit 
simulation considering the presence of moons and planetary rings around the planet. 
The moon's orbit is considered to be circular and coplanar with the planetary orbit. 
The other physical and orbital parameters of the star, planet, moon and rings can be 
adjusted in each simulation. It is possible to simulate as many successive transits as 
desired. Since the presence of spots on the surface of the star may produce a signal 
similar to that of the presence of a moon, our model also allows for the inclusion of 
starspots. The result of the simulation is the light curve with the planetary transit. 
White noise may also be added to the light curves, to produce curves similar to those 
obtained by the CoRoT and Kepler space telescopes. The goal is to determine the 
criteria for detectability of moons and/or ring systems using photometry. The results 
show that it is possible to detect moons with radii as little as 1.3 R© with CoRoT and 
0.3 R® with Kepler. 



Subject headings: extrasolar planets, exomoons, planetary transits 
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Introduction 



Since the detection o f the first extra-solar planet orbiting a solar- like star in 1995, 51 Peg 
b (jMavor &: Queloa Il995l ). several new discoveries have been made. The most used planetary 
detection methods are radial velocity and transit. Together, both techniques led to the discovery 
of the vast majority of the over 570 planets observed so far (according to the Extra-solar Planets 
Encyclopaedia, |http://ex oplanet.eu). When the same planet is detected by both methods, the 
complete set of planetary and orbital parameters is determined. Two space missions designed to 
detect planetary transits are operational right now: CoRoT, from the CNES with a consortium of 
many countries, and Kepler, from NASA. 

One of the motivations for the search of extra-solar planets is the quest for places suitable for 
the emergence of lif e. For this to be p o ssible, the planet ne eds to be located in a place known as 
the habitable zone (jHuang. S. S.I Il959i ; iKasting et al.l Il993l ) . This is the region around the host 
star, where the radiation received by the planet is such that allows for water to exist in its liquid 
state. Furthermore, the planet needs to be rocky. Due to the observations bias, the majority of the 
planets found up to date are gaseous giants, what excludes them from being habitable even if they 
are in the habitable stellar zone. 

Through observation of the Solar System, we can infer that the existence of moons and plan- 
etary rings around gas exoplanets may be common. Also, we can deduce that moons are rocky by 
nature. Planetary rings may exist in very different configurations, from the rarefied rings of Uranus 
to the giant rings of Saturn. The detection of both moons and/or rings may uncover important 
clues about the planet formation and migration process. Moreover, if a gas planet located in the 
stellar habitable zone has moons, these moons may also be habitable, making them a good place 
for life to evolve. 



Sartoretti k, Schneider! (119991 ) suggest that exomoons could be detected using photometry. 
The presence of moons around an exoplanet would cause distortions in the transit light curve, in 
the form of steps or asymmetries. The detection of these distortions would provide the moon's 
orbital parameters. Using this technique and photometric data from the Hubble Space Telescope, 
Brown et al.1 (|200ll ) s earched for mo ons and rings around the exoplanet HD 209458b. However, 



nothing was found. iKippineJ (|201ll ) proposes an algorithm capable of generating transit light 



curves of a planet with a moon analytically, and shows how th is can be used on the search for 
moons using photometry. On the other h and. lHan fc Hani (120021) st udied the detection possibility 
using gravitational micro- lensing events. Williams &: Knackel ( 2004 ) showed that it is possible to 
detect moons in orbit of gaseo us giant planets in the habitable zone using spectroscopy in the 1.5 
- 5 /xm band. iKipping d2009al lbh suggests that moon detection is possible from the effects on the 
planet motion, that may be detectable through precise measurements of the timing and duration 
of the transits. These are known as transit time variation (TTV) and transit duration variation 
(TDV). 



Barnes fc Fortnevl (|2004l ) suggest that rings around a Saturn-like exoplanet could be detected 



- 3 - 



with photometric precision of (1 - 3) x 1CP 4 and time resolution of 15 minutes, as long as the 
ring is not viewed close to edge-on. This resolutio n is within t he ac curacy that the Kepler Space 
Telescope can achieve for Sun-like or brighter stars. lOhta et al.l (|2009l ) study ring detection through 
photometry and spectroscopy, and conclude that ri ngs with significantly great inclination can be 
detected by today radial velocity precision of 1 m/s. ISchlichting fc Chang! (j201ll ) study the nature 
of rings that can exist around exoplanets, and conclude that rings around close-in planets can show 
non-trivial Laplacian planes. This is caused by the increased effects of the planet's quadrupole 
moment. These so called warped rings furnish important clues about the planet interior and spin 
period. 

Despite all these proposed methods, no exo-moon or planetary rings have been observed so 

far. 

Here, we propose a model of planetary transit simulation that may be used in the search for 
exo-moons through their signature in the light curve of their host planets. Beside the moons, the 
model may also simulate the transit of a planet with a ring aro und it. T h e mod el we propose here 
is capable of generating transit light curves numerically, unlike iKippina (|201ll ). Albeit the larger 
processing time needed to fit light curves using this method, it has the advantage of being easily 
adaptable to include new features, as additional moons, changes in the ring plane or starspots. The 
goal is to apply this model to the observations of the CoRoT and Kepler space telescopes, looking 
for the photometric signals that may indicate the presence of moons or rings around exoplanets. 



Our model does not consider the detection of the warped rings shown by ISchlichting &: Chang 



(|201ll ). Such rings are predicted around close-in exoplanets, with orbital period of a few days. 
Thus, our model is limited to the detection of rings located in the planet's equatorial plane, with 
larger semi-major axis. As the occurance of such warped rings depends on the quadrupole moment 
J2 and on the density of the planet, it is difficult to stablish a general criterion to the validity of our 
model. For example, for the case of an exoplanet with densit y 3 g cm" 3 and J? = 1 ~ 2 , w arped rings 



occur for semi-major axis below 0.05 AU (see Figure 6 on ISchlichting &: Chang) (120111 )). Planets 
with larger semi-major axis will have rings located in the planet equatorial plane and could be 
detected by our model. 

The next section describes the model used in this work, whereas Section 3 presents the appli- 
cation of this model to a few study cases. Section 4 discusses the detectability threshold of moon 
and ring systems. Finally, the last section lists our main conclusions. 



The model 



The model used here is based on that of ISilval (|2003l ). where the star is considered a disc with 
limb darkening and the planet a completely dark disc. Both stellar and planetary parameters can 
be fit by the model. To minimize the number of free parameters, we consider a simple model where 
the moon's orbit is coplanar with the planetary one, and is also circular. Looking at the Solar 
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System's moons, we see that these assumptions are a little limiting. However, as this is our first 
study on the detection of moons, we choose to consider the simplest model. In a future work we 
can consider a more realistic model, including orbital inclinations and eccentricity fo r the moon. 



Regarding to the planet eccentricity, following the strategy proposed by iKippinel (|2008l ). the model 
simulates eliptical orbits only for the cases were eccentricity was measured by other methods, such 
as radial velocities. 

Stellar input parameters to the model are: rotation period, radius, and mass (the latter two in 
solar units). The rotation period is used in the calculation of the longitude position of starspots in 
successive transits. We also have two limb darkening parameters. The planetary input parameters 
are: orbital period, radius, mass, orbital inclination angle, eccentricity and position of the peri- 
astron. From the orbital period, the orbit's semi-major axis is calculated by Kepler's Third Law. 
Lastly, the moon input parameters are orbital period, radius, mass, and the angular position of 
the moon in its orbit at the starting time of the first transit. For the rings, the input parameters 
consists of inner and outer radii, two inclination angles, one with respect to the plane of the sky 
and the other with the planet's orbital plane. The transparency of the rings may also be fit. 

The temporal resolution of the resulting light curve is chosen within the model, by adjusting 
the time interval for each simulated photometric data point. Here we used a temporal resolution 
of 32 s for future comparison with CoRoT data, and 1 min for Kepler data. 

Using stellar and planetary parameters, the planetary orbit is calculated. The line of transit is 
the projection of this orbit on the stellar surface. At each time interval, the program calculates the 
position of the planet and its moon. The total luminosity of the star-planet-moon system is then 
determined by summing the intensities of all the pixels in the image, where the pixels occupied 
by the planet and its moon have zero value. In the case of rings, the occupied pixels intensity is 
multiplied by the ring transparency factor. At the end of the transit, the light curve is normalized 
to one from the outside transit data points. 



Figure Q] shows two examples of transits from the model. Figure 1(a) shows the transit of a 



planet and a moon in front of the star, and Figure 1(b) shows the transit of a planet with rings 



It is also possible to add spots on the stellar surface, where each spot is modeled by three extra 
parameters: radius, intensity (with respect to stellar central intensity), and longitude position at 
the start of the first transit. The spot latitude is considered to be that of the projected transit line. 
From the stellar rotation period, the spot longitude position on consecutive transits is calculated. 
A further assumption is that the spot position does not change within a single transit (usually 
lasting a few hours), which is reasonable when the stellar rotation period is much larger than the 
transit duration. 
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Fig. 1. — (a) A Jupiter-like planet and a 2 R^ moon transiting in front of the star, (b) A Saturn- 
like planet with rings transiting in front of the star. The rings have the same dimensions as Saturn's. 
The angle with the sky plane is 78°, and the angle with orbital plane is 20°. Both images were 
obtained from the model. 



2.1. Planet-Moon dynamics 

In this work, no attempt has been made for a dynamical analysis of a plane t with moons or 



rings 



We used for this matter, regarding to planet-moon systems, the results of iDomingos et al 



(|2006l ). where the authors through numerical simulations obtain semi empirical expressions for the 
semi- major axis maximum and minimum that a satellite in orbit of a giant e xtras olar planet may 
have as a function of the stellar and planetary parameters. IDomingos et al. also estimates 



the maximum mass of a satellite such that its orbit is stable. An interesting result is that planets 
with very close in orbits to their host star can only hold moons with very small masses, but more 
distant planets may have stable moons with masses even larger than that of the Earth. 

For example, Rho Crb at a distance of 0.224 AU from its host star, with an orbital period 
of 39.8 d, may harbor a satellite of up to 1.5 M®. Further away planets, such as HD 92788, at 
approximately 1 AU from its star, allows for a satellite with a maximum mass of 1 x 10 4 M®. 
The transiting extrasolar planets with the longest periods observed so far are CoRoT-9b (95.27 d), 
HD 80606 (114.34 d) and Kepler-11 g (118.37 d). Thus, these planets may harbor a moon with 
mass similar or even larger than Earth, that may be detectable from the application of the present 
method to CoRoT or Kepler lightcurves. 
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2.2. Planetary transit 



Figure [2] shows two examples of lightcurves ob tained from our m odel. The figure on the left, 
2(a)| is for the transit of a planet like HD 20945 8b (Southworth ( 20101 )). whereas the figure on the 
right, |2(b)[ is for a planet similar to CoRoT-2b (|Alonso et al.1 (|200sl )). 



(a) 



(b) 



Fig. 2.— Transit model of two planets, (a) HD 209458b (R = 1.38 Rj, P = 3.524 days, i 
and (b) CoRoT-2b (R = 1.465 Rj, P = 1.743 days, i = 87.84°). 



86.59° 



2.3. Transit of a planet with a moon 



Figure [3] shows two examples of lightcurves of a planet with an orbiting moon. Both panels are 
for a Jupiter size planet transiting a star identical to the Sun with an orbital period of 45 days. In 



the model, the added moon around this planet has 2 Earth radii. In 3(a) , the moon orbital period 
is 5 days, while in Figure |3(b)| it is 2 days. Also, the position angle of the moon with respect to 
the planet is 0° in the first case, and 60° in the second one. 




(a) 



(1)) 



Fig. 3. — Model light curve of a transiting hot Jupiter (R = 1 Rj, P = 45 days, i = 89.5°) 
around a Sun-like star, and an orbiting moon with 2 R$. The moon period and position angle are, 
respectively: (a) Pm = 5 days and ip = 0°; (b) Pm = 2 days and (p = 60°. 
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When adding a moon around a planet in the model, two characteristics signatures may occur. 
The first one is a step, that may appear at the beginning or the at the end of a transit, while the 
second effect is a discontinuity, a jump, seen in the bottom of the transit. The step at the beginning 



(see Figure 3(a) ) occurs when the moon ingress happens before that of the planet, while the step 
near the end of the transit (Figure [3 (b)[ ) is caused by the moon egress after the planet. Both effects 
are marked with an upward arrow. 

Furthermore, when the moon orbital period is small, an asymmetry arises in the bottom of the 
transit, shown by the downward arrows in Figure [3l This asymmetry is generated when the moon is 
eclipsed by the planet. It may also happen if at the start of the transit the moon is already behind 
(or in front) of the planet, but falls out of eclipse during the transit. The width and intensity of 
this signal depends on the moon's relative orbital period, radius, and also on the position angle ip 
of the moon with respect to the planet at the initial transit time. 



2.4. Transit of a planet with rings 

Figure H] presents the modeled light curve produced by a transiting Saturn-like planet around 
a star with solar parameters (left panel). To enhance the signal produced by the ring presence, 
their size were augmented by 50% (right panel). 

The main effect of the planetary ring system is to increase the covered area of the stellar 
surface. The result of which is to deepen the transit light curve. Moreover, the transit shape 
becomes rounder at the ingress and egress instants, than if no rings are present. Both these factors 
depend strongly on the ring transparency, being larger in the case of more opaque rings, as expected. 

An interesting effect appears when the ring size is increased. In this case, the light curve 
becomes more triangular, similar to the light curves produced during grazing eclipses of binary star 
systems. 



Ti™ (I,) Tim, (h) 

(a) (b) 

Fig. 4. — Transiting Saturn-like planet (R = 0.84 Rj, P = 5 days, i = 88°) in front of a solar-like 
star. Planetary rings with sizes (a) equal and (b) increased by 50% with respect to those of Saturn. 
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2.5. Starspots 



If the star has one or more spots on its surface, similar to sunspots, ISilval (|2003l ) has shown 
that the photometric light curve during a planetary transit will display a small "bump". This 
occurs because when the planet occults a darker region of the stellar surface, occupied by the spot, 
there is a relative increase in the total intensity Moreover, the detection of the same s t arspo ts on 
different transits may be used to determine the rotation period of the star (jSilva- Valid (120081 ) ). 

Figure [5] shows the same light curves from Figure [21 except now a spot was added to the stellar 
surface. This spot, located at 30° longitude (measured from the center of the star), has a size of 
half each planet radius and a relative intensity of 0.6 times the stellar central intensity. 



(a) 



(b) 



Fig. 5. — Planetary transit light curves for (a) HD 209458b and (b) CoRoT-2b, where the "bump" 
in the bottom of the transit indicates the presence of a spot on the stellar surface. In each case, 
the spot has half the planet radius and intensity 0.6, and is located at 30° with respect to the star 
center. 



In some cases, the photometric signal produced by the presence of moons may be mistaken by 
that caused by starspots. Distinguishing between the two cases is possible when successive transit 
are observed, since the spots signal will have a periodicity equal to the stellar rotation period, while 
those of the moons will have a different period, following its own orbital period. 

Examples of how confusing things can get, especially how the presence of starspots may hinder 
the detection of moons, are shown in Figure El The figure presents the transit of a planet with an 
orbiting moon in front of a star, in which surface there is a spot, for two situations: only one spot 
(6(a) ) and star with 4 spots ( |6(b)[ ). In both panels, the planet is identical to Jupiter with a period 
of 45 d and the moon has 3 Earth radius and an orbital period of 2 d around the planet. In the 
case of Figure |6(b)[ it is clearly seen that the presence of starspots can significantly difficult the 
detection of any moon effect on the light curve. 
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(a) (b) 

Fig. 6. — Jupiter-like transiting planet (R = 1 Rj, P = 45 days, i = 89.5°) with a 3 R^ moon around 
a star equal to the Sun. (a) The star has a spot with 0.5 planet radius, intensity 0.6, located at 
30°. (b) The star has 4 spots, with radius respectively 0.7, 0.6, 0.7 and 0.8 planet radius, intensity 
0.6, 0.5, 0.6 and 0.7, located at -30, -15, 5 and 30°. Each spot signal appears duplicated, once when 
the moon occults it and the other when the planet eclipses the spot. 



3. Test results from the model 



To verify the detectability of moon and ring systems around a planet, it is possible to add 
Gaussian noise to the modeled light curve and then try and recover the input parameters from 
a fit. For CoRoT data simulation, we used a noise level of 6 x 10 -4 , such as that observed for 
CoRoT-2 (jSilva-Valio k. Lanzall2010l ). and in the case of the Kepler data, 2 x 10~ 5 . 



Once the noise has been added to the mod el light curve, we fit t his light curve with the following 
sequence of algorithms. First the PIKAIA dCharbonneaul 1 19951) algorith m is used to obtain a 
first guess that is then inputed into the AMOEBA (jNelder k. Meadlll965l ) routine. PIKAIA is a 
genetic algorithm the main advantage of which is to find the global maximum of the function to 
be maximized, using only the allowed interval value for each free parameter. These results are 
then used as input to the AMOEBA routine that refines the result using the simplex method for 
convergence minim ization. The erro r estimates for the fitted parameters were obtained using the 
MPFIT algorithm (|Markwardtll2009h . 



After the fitting procedure, we compute the goodness- of -fit using the parameter Q, as defined 
in IPress et al.1 ( 1992 ) as the probability that the chi- square should exceed a particular value \ 2 by 
chance and calculated by: 

2' 2 



Q(x 2 W) 



where P is the incomplete gamma function P(a,x). Generally, if Q is larger than 0.01, then the 
fit is acceptable. In this case, it can be said that such event (moon or ring) is detectable. If Q is 
between 0.001 and 0.01, then the fit may or not be acceptable, and more studies about the noise 
are required. This can be caused, for example, if the measurement errors are underestimated. Q 
below 0.001 indicates a bad fit. On the other hand, Q very close to 1 is not expected also. This is, 
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literally, too good to be true. This can be caused for nonnormal error distribution, or overestimated 
measurement errors. Next, we show some examples of the fit to light curves with noise added, and 
the results they yield. 



3.1. Fit of a planetary transit 



The first test to the model is to apply it to the two transiting planets described in Section \2.2\ 
where Gaussian noise has been added similar to that of CoRoT and Kepler observations. In this 
fit, all the stellar parameters and the orbital period (or semi-major axis) of the planet were fixed, 
and the only free parameters were the planetary radius and the orbital inclination angle. 

The comparison between the original parameters and those obtained from the fit of the noisy 
modeled light curve are seen in Tables [1] and [2j and and in Figures [7(a)] and [7(b)] for HD 209458 
and CoRoT-2. Both light curves had noise level similar to the CoRoT satellite. The two panels 
show the simulated data points as dots and the fit as a solid line. The residuals are shown in the 
lower part of each panel. 





(a) 



(b) 



Fig. 7. — Light curve fit of model planetary transit for (a) HD 209458b and (b) CoRoT-2b. The 
dots represent the noisy light curve data points, whereas the solid line shows the best yield fit. 



Table 1: Fit results for planetary transits of HD 209458b and CoRoT-2b for CoRoT simulations. 



Planet 


Parameter 


Original Value 


Fit Result 


Reduced \ 2 


Q 


HD 209458b 


Planet radius (Rj) 
Orbital inclination (°) 


1.38 

86.59 


1.382 ± 0.025 
86.6 ± 0.3 


1.0639 


0.098 


CoRoT-2b 


Planet radius (Rj) 
Orbital inclination (°) 


1.465 
87.84 


1.466 ± 0.016 
87.84 ± 0.21 


1.0157 


0.3890 
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Table 2: Fit results for planetary transits of HD 209458b and CoRoT-2b for Kepler simulations. 



Planet 


Parameter 


Original Value 


Fit Result 


Reduced x 2 


Q 


HD 209458b 


Planet radius (Rj) 
Orbital inclination (°) 


1.38 
86.59 


1.3797 ± 0.0011 
86.594 ± 0.015 


1.0627 


0.1747 


CoRoT-2b 


Planet radius (R,/) 
Orbital inclination (°) 


1.465 
87.84 


1.4636 ± 0.0007 
87.839 ± 0.009 


1.1156 


0.0859 



3.2. Transit fit of a planet with moon 

If the planet has a moon in orbit around it, there are two factors that make its detection 
difficult. The first one is the noise, as expected. In cases where the moon is very small, their 
signature in the light curve is flooded by the noise level, leaving it undetectable. Another difficulty 
is the presence of starspots on the stellar surface, mentioned above. 

To demonstrate how moons can be detected, we consider the case of a hot Jupiter with a 
period of 45 days around its host star, modelled with no spots. In orbit around this planet there 
is a moon with 3 Earth radii and a period of 5 days. Simulations were performed for both noise 
levels of CoRoT and Kepler. The parameters obtained from the fit and the original ones are listed 
in Table [31 whereas the fit results are plotted in Figure [U for CoRoT noise level. The fit was 
performed for 5 consecutive transits, but the figure shows only the first. 



Fig. 8. — Light curve fit for the first transit of a hot Jupiter (R = 1 Rj, P = 45 days, i = 89.5°) 
with a moon (Ry[ = 3 R®, Py^ = 5 days) considering the noise level of CoRoT. 



We noticed that the larger the number of transits considered in the fit, the smaller the number 
of iterations necessary for the fitting algorithms to obtain a good fit, that is, the easier it is to 
detect the moon. 
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Table 3: Fit result of 5 transits of a planet with an orbiting moon. 



Noise Level 


Parameter 


Original Value 


Fit Result 


Reduced x 2 


Q 




Planet radius (Rj) 


1.0 


0.9976 ± 0.0005 






CoRoT 


Moon period (days) 


5.0 


5.001 ± 0.012 


1.0070 


0.3954 








Moon radius (R©) 


3.0 


2.979 ± 0.019 








Moon position angle (°) 


0.0 


0.21 ± 0.11 






Kepler 


Planet radius (Rj) 
Moon period (days) 
Moon radius (R©) 
Moon position angle (°) 


1.0 
5.0 
3.0 
0.0 


0.9999 ± 0.0007 
5.000 ± 0.023 
2.999 ± 0.024 
-0.4 ± 0.6 


1.0419 


0.1337 



3.3. Transit fit of a planet with rings 

As seen above, planetary rings can be detected by the round shape of the transit at ingress 
and egress times, besides the larger depth of the transit (due to the increase in the stellar surface 
area covered). For these signatures to be detectable, the noise is the predominant factor. If the 
signature is too small, the noise will completely overcome it. 

To exemplify, we simulate the transit of a planet identical to Saturn and its rings, with 30-days 
period and 89° orbital inclination. The transit is modeled for two cases: that with the noise level of 
CoRoT and the other with that of Kepler telescope. In both cases, the light curve was submitted 
to a fit process assuming only the presence of the planet. The result of this preliminary fit is shown 
in Figure [9] and Table [U 

Table 4: Fit results for a Saturn-like planet simulation, considering only the planet. 



Noise level 


Parameter 


Original Value 


Fit result 


Reduced x 2 


Q 


CoRoT 


Planet radius (Rj) 
Orbital inclination (°) 


0.84 
89 


0.934 ± 0.022 
89.0 ± 0.4 


1.0291 


0.2050 


Kepler 


Planet radius (Rj) 
Orbital inclination (°) 


0.84 
89 


0.9434 ± 0.0009 
88.999 ± 0.019 


25.10 






From the analysis of the residual plotted on Figure [9] and the Q value in Table dj we see 
that the fit considering only the planet yields a satisfactory result only for the CoRoT noise level 
(Figure [9(b)[ ). This indicates that a ring system such as that of Saturn could only be detected by 
the Kepler satellite. The CoRoT satellite would identify this light curve as that produced by a 
giant planet with a radius slightly larger than that of Saturn (12.3 %), but without rings. 

Next, we fit only the light curve plotted in Figure [9(b)] but now considering the presence of 
rings, that is the ring parameters are taken into account. The resulting fit parameters are listed in 
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(a) 



| 0.0002 | 
.1 0.0000 § 
K -0.0004 | 




-6 -4 



Time (h) 

(b) 



Fig. 9. — Saturn-like planet (R = 0.84 Rj, P = 5 days, i = 88°) transit fitted without considering 
the presence of its rings for noise levels of the (a) CoRoT and(b) Kepler telescopes, respectively. 



Table [5l whereas the result of the fit is plotted in Figure [TOj Comparing the residual of Figures |9(b)| 
and llOl one can see that the fit is much better for the model where the rings were included. 





\ j 















Fig. 10. — Saturn-like planet (R = 0.84 Rj, P = 5 days, i = 88°) transit fitted, considering now 
the presence of rings. 



The Q value in Table [5] and the residuals in Figure [TU] show that the fit was satisfactory. 



4. Detection limits 

To determine the radius of the smallest detectable moon, several simulations were performed 
decreasing the radius at each run. Every simulation had its light curve fitted yielding the moon 
parameters. Whenever the Q value indicates a good fit, we consider the moon detectable. Thus, in 
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Table 5: Resulting parameters for the case of a transiting Saturn, as would be observed by Kepler. 



Parameter 


Original Value 


Fit Result 


Reduced \ 2 


Q 


Planet radius (Rj) 


0.84 


0.839 ± 0.005 






Orbital inclination (°) 


89 


89.00 ± 0.03 






Inner radius (Rp) 


1.11 


1.11 ± 0.06 










1.0173 


0.3553 


Outer radius (Rp) 


2.32 


2.316 ± 0.026 






Angle with the sky plane (°) 


78.0 


77.9 ± 1.5 






Angle with the orbital plane (°) 


20.0 


19.9 ± 0.5 






Transparency 


0.5 


0.51 ± 0.06 







the case of CoRoT, the smallest detectable moon in orbit around a Jupiter size planet has a radius 
of 1.3 R®, while in the case of Kepler, this minimum radius decreases to 0.3 R®. The latter moon 
would be smaller than Ganimedes, the larger moon of our Solar System, which size is 0.41 R®. In 
all the simulations, the star was considered Sundike with no spots. 

For the case of planetary rings, the process is slightly different, because, in this case, the 
biggest problem is distinguishing between single planets and those with rings. In this case, then, 
we adopt a criterion that takes into consideration the residuals. When fitting the light curve of a 
planet with rings, first only the planet parameters are fit. If the Q value indicates a bad fit, then 
we may conclude that the rings can be detectable. 

To establish a detectability criterion, first we define the effective area of the rings as the ring 
projection on the plane of the sky multiplied by its opacity, according to: 

A eff = (R 2 E -Rj) it cos i/j(1-t) (2) 

where Ri and Re are the internal and external radii, respectively, ip is the inclination angle of the 
rings relative to the plane of the sky, and r is the transparency. Totally transparent rings have 
r = 1, whereas totally opaque ones have a zero value of transparency. 

The results of our simulations showed that in the case of the Kepler data, planetary rings are 
detectable when the effective area corresponds to approximately 3% of the area covered by the 
planet. For example, rings with r = 0.5 and inclination angle of 78° with the plane of the sky in 
orbit a Saturndike planet, will be detectable with rings about 51% the size of those of Saturn. For 
the CoRoT data, the effective area corresponds to approximately the area of the planet. Therefore, 
for the same planet as in the previous case, the rings should have sizes 45% larger than those of 
Saturn to be detectable above the CoRoT noise level. 
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5. Conclusion 



Here we reported the results of a model capable of simulating planetary transits, including the 
effect of existing moons and rings around that planet. Light curves were presented for each case. 
Estimates of the threshold detection for such events were made. It was shown how these events 
may be detected using the PIKAIA and AMOEBA routines for curve fitting. As expected, the 
main limiting factor for the possible detection is the noise level of the telescopes. This can be seen 
by comparing the results from simulations of photometric light curves of CoRoT and Kepler. 

The CoRoT data allows for the detection of moons with at least 1.3 R® in size around a planet 
identical to Jupiter. In the case of Kepler data, this detection limit on the moon radius drops to 
0.3 R®. For planetary rings, a Saturn-like planet would be detectable with Kepler, however, this 
same planet-ring system would be mistakenly inferred as a single planet with slightly larger radius 
by CoRoT. Detection limits for both telescopes were estimated above. 

Next, this model will be applied to real light curves to search for the types of signatures 
described in this work. Data from both CoRoT and Kepler will be used. An improvement that 
can be made to the model is the implementation of analysis of the transit time and duration 
variat ions caused by the presence of the moon that makes the planet orbit the common center-of- 



mass 



kippinglboOSy lbh. These timing variations would confirm the presence of moons. Finally, by 



detecting the moon by both methods, photometry and timing variation, it is possible to establish 
a complete set of physical parameters for the satellite. 
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